Temperature determined by isobaric yield ratio in heavy-ion collisions 
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Background: Temperature (T) in heavy-ion collision is an important parameter. Previously, many works have 
focused on the temperature of the hot emitting source. But there are few systematic studies of the temperature 
among heavy fragments in peripheral collisions with incident energies near the Fermi energy to a few A GeV, 
though it is very important to study the property of neutron-rich nucleus in heavy-ion collisions. 

Purpose: This work focuses on the study of temperature associated with the final heavy fragments in reactions 
induced by both the neutron-proton symmetric and the neutron-rich projectiles, and with incident energy ranges 
from mA MeV to lA GeV. 

Methods: Isobaric yield ratio (lYR) is used to determine the temperature of heavy fragments. Cross sections of 
measured fragment in reactions are analyzed, and a modified statistical abrasion-ablation (SAA) model is used 
to calculate the yield of fragment in 140A MeV ^^Ki + ''Be and lA GeV ^^®Xe + ^°*Pb reactions. 

Results: Relatively low T of heavy fragments are obtained in different reactions {T ranges from 1 to 3MeV). 
T is also found to depend on the neutron-richness of the projectile. The incident energy affects T very little. 
Afi/T (the ratio of the difference between the chemical potential of neutron and proton to temperature) is found 
to increase linearly as N/Z of projectile increases. It is found that T of the **Ca reaction, for which lYRs are 
of A < 50 isobars, is affected greatly by the temperature-corrected AB{T). But T of reactions using lYRs of 
heavier fragments are only slightly affected by the temperature-corrected AB(T). The SAA model analysis gives 
a consistent overview of the results extracted in this work. 

Conclusions: T from lYR, which is for secondary fragment, is different from that of the hot emitting source. T 
and A/i are essentially governed by the sequential decay process. 

PACS numbers: 25.70.Pq, 21.65. Cd, 25.70.Mn 



I. INTRODUCTION 23 

24 

In heavy-ion collisions (HIC) above the intermediate 
energy (with incident energy > 20 A MeV) , the temper- ^'^ 
ature (T) is high enough to provide the environment for " 
nuclear liquid-gas transition. Many works have focused 
on the critical point in T experimentally [ll-Q and the- 
oretically [1, [l^l . T is also important in determining 
the symmetry energy of a neutron-rich nucleus pdl - fTHj . 
Due to the complexity of the HIC processes and the im- 
portance of T, many methods have been developed to 
obtain T from thermal energy p^. excitation energy 
[20j . isotopic yields (Al bergo t hermometer) [sl. [2T| - [26| . and 
kinetic energy spectra [27H29| . T based on these methods 
have differences but can be related to each other [l^, [2^ . 

The yield of fragment in HIC, on the one hand is 
greatly influenced by T, and on the other hand, is deter- 39 
mined by its free energy due to non-zero T [l[ . Thus the 40 
yield of fragment can constrain both its binding energy 41 
and T. In some models which estimate isotopic yield, T 42 
is an important parameter (30l-l33| , while in other works 43 
the isotopic yield is also used to estimate the binding en- 44 
ergy of unknown isotopes jsof . In works using the yield of 45 
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fragment to constrain the binding energy, for example, in 
the study of the symmetry energy of fragments, it is pro- 
posed that the isobaric yield ratio cancels out the energy 
term which only depends on the mass of fragment, thus 
the specific energy term can be extracted [SJ, HH . But 
the shortcoming of this method is that the coefficient of 
energy-term and T can not be separated, thus they must 
be viewed as whole parameters, for example, in the study 
of the symmetry-energy coefficients to T (asym/T) of the 
neutron-rich nucleus [sj, [36^38} . If the binding energy of 
fragment at nonzero T is known, the yield of fragment 
can also be used to determine T. One point to remember 
is that the measured heavy fragment in HIC undergoes 
the sequential decay and deexcitation processes, which 
makes it cool down. T of a heavy fragment should be 
quite lower than that of the hot emitting source [13, ■ 
Based on the free energy, the modified Fisher model 
(MFM) well describes the isotopic yield distributions of 
intermediate mass fragments produced in proton-induced 
multifragmentation at relativistic energies P, H^. The 
MFM has been used to study the behavior of fragments 
near the critical point of the liquid-gas transition 
Q. In the MFM, the yield of an isotope is determined 
by chemical potential, free energy and entropy. In this 
article, we will use the isobaric yield ratio to constrain 
T in HICs. The article is organized as follows: First, 
we briefiy introduce the isobaric ratio method based on 
MFM. Second, we verify that the binding energy of nu- 
cleus at zero temperature can be used instead in the iso- 
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baric ratio method. At last, T determined by isobaric 45 
yield will be shown and the results will be discussed. « 



II. ISOBARIC YIELD RATIO METHOD 

50 

Following the MFM [ll, li^l , the generalized expression 51 
of the yield of a fragment with mass A and neutron-excess 
I{I = N - Z) is 

Y{A, I) = CA-^exp{[W{A, I) + ^„iV + fipZ]/T 

+Nln{N/A) + Zl7i{Z/A)}, (1) 

where C is a constant. The A^"^ term originates from 
the entropy of the fragment, r's for all fragments are 
identical. /i„ and are the neutron and proton chemi- 
cal potentials, respectively, and /) is the Helmholtz 
free energy of the cluster (fragment). In principle Eq. ([1]) 
should be applied to hot nuclear matter near the critical 
point. However, when Eq. ([T]) was applied to the cold 
fragments in Ref. [4l[, they showed that useful infor- 
mation can be extracted to elucidate the effects of the 
secondary decay process. In that case, T and other pa- 
rameters do not correspond to those in the primary hot 
nuclear matter, but do correspond to those modified by 
the secondary decay process. Therefore it is still use- 
ful to apply Eq. ([T]) to the experimentally observed cold 
fragments to elucidate the effect of the sequential process 
on the characteristic physical parameters, such as T or 
chemical potential, as discussed below. 

Defining the yield ratio between isobars differing by 2 
units in /, we have 

Ril + 2, /, A) = Y{A, I + 2)/Y{A, I) 
= exp{[W{I + 2,A)-W{I,A) + A^t]/r 

+ s„,ui + '2,A)~s„,,M,A)}, (2) 

where Smi4l,A) = Nln{N/A) + Z\n{Z/A), and A/z = 
— Mp- Taking the logarithm of R{I + 2, 7, A), one 
obtains 



lni?(/ + 2, /, A) - AS* = (Aiy + Afi)/T 



(3) 



where AS" = ^mixU + 2,A) - S'mix(/,A), and AW =i 
W(J + 2,A)- W{I, A) is the difference between the free 1 
energies of isobars. W{I,A) is supposed to equal the 1 
binding energy (B) of a fragment at given T and density < 
p [3 [written as B{p,T)]. If B{p,T) (which includes the < 
contributions from the binding energy and entropy) is t 
known, T and A/i can be obtained using the isobaric t 
yield ratio (lYR) from Eq. ©. < 
The known B of the nucleus is for T = [written as i 
B{Q)]. When T 7^ 0, entropy contributes to the binding t 
energy [written as B{T)] and makes B{0) ^ B{T). The. 
other good news is that since AW can serve as the in- < 
dependent variable, actually B{T) is not important any- 1 
more in Eq. ([3]). If B{0) can be used instead of B(T) in 1 
the calculation of AW, the temperature of the fragment 1 
can be extracted using Eq. (jS]). Thus before determin- 1 
ing T using the lYR method, whether AB{0) = AB{T) ■ 



should be evaluated [AB{T) = S(/+2)(T) - S(/)(T) is 
the difference between the binding energies of the 1 + 2 
and / isobars]. 

Using the density-functional theory based on the 
Skyrme interaction (SKM), the T dependence of the 
binding energy of a finite nucleus has been proposed to 
be lii, 



B{A,I,T) 



(15.31 - 0.04r2)^ 
(18.30 + 0.50T2)A2/3 
(19.69 + 0.42T2)/2/^ 
- (33.18 + 2.06T2)/2/A4/3 + 



^2 



^4/3 



+ A(iV,Z), (4) 



where Edi / and E^x are the coefficients for the diffuseness 
correction and the exchange correlation to the Coulomb 
energy. A(iV, Z) is the pairing-energy term. 



AB=B,„„(T)-B,„(T) 
• T=1 ▲ T=2 ♦ T=3 
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FIG. 1. (Color online) Ato = [AB(T) - AB(0)]/AB(0) of 
isobars with I from -1 to 11 at different T. AB{T) is the 
difference between the binding energies of the 1 + 2 and / 
isobars at T. AB{T) = B(/+2)(r) - B^^i){T). 



The binding energy of the nucleus at T from to 5MeV 
are calculated using Eq. (|4]). To see how fast AB{T) be- 
tween isobars increases with temperature, the values of 
Ato = [AB{T) - AB(0)]/AB(0) arc plotted in Fig. [1] 
At T ~ IMeV, Ato of the /= -1 - 11 isobars are 
very close to zero. Ato decreases as the mass becomes 
larger, i.e., the larger the A of nucleus, the closer Ato 
is to 0. For isobars of small mass, Ato decreases simi- 
larly as T decreases. At T > 3MeV, the A < 50 isobars 
show relatively large Ato (Ato < —0.1 corresponds to 
uncertainty larger than 10%). AB{T) of A < 50 isobars 
should be used more carefully around T > 3McV. In Eq. 
(jSI AW is the difference between free energies of isobars, 
thus the very little Ato for I > 5 and A > 50 isobars 
occurs when AB{T) is replaced by Ai3(0), and we need 
not know the actual B{T). The smaller Ato of isobars at 
T, the closer T obtained from Eq. ^ approximates the 
real value. Thus theoretically, the yield ratios of / > 5 
isobars are suitable obscrvables to extract T. Replacing 



3 



the AW term in Eq. ([3]) by AB{0) (hereafter denoted as 
AB), one obtains, 



a: 



lni?(/ + 2, 1, A) - AS = {AB + A^i) /T 



(5) 
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FIG. 2. (Color online) The correlation between lYR and AB 
of isobars in the 6iA MeV *^Kr + ^Be projectile fragmenta- 
tion reaction [ij. The lines are the fitting results of lYRs 
using Eq. ([S}. The solid line represents the fitting line of all 
the data. The inserted figure shows the fitted values of T, in 
which the solid line represents T fitted from all the data, and 
the dotted line represents T =2.2 MeV used in Ref. [s^]. : 
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FIG. 3. (Color online) The correlation between lYR and AB 
of isobars in the following reactions: (a) \A GeV ^"^"^Xe + Pb 
[H^b) \A GeV i^^Xe -f Pb [H, (c) 790A MeV '^''Xe + 
Al m, (d) \A GeV "'Sn+^'Sn and i24sn+i24gj^ g^j^ ^g-, 
140A MeV **Ca + ^Be/^^^Ta [l!], and (f) f404 MeV '^''Ni 
+ "Be/^^^Ta [H]. AB is calculated using the experimental 
binding energy of the nucleus in Ref. [ij • The lines denote 
the fitting results using Eq. ((SJ and T is the temperature 
obtained. 



T = 2.23MeV. By analyzing the results of T, conclusions 
below can be drawn: 



III. RESULTS AND DISCUSSIONS 

By analyzing lYR in HIC, T and A/i can be determined 
using Eq. ([5]). The analysis is performed using the lYR 
of / > 5 isobars. First, the dependence of T on / is 
investigated using the data in the 64^ MeV ®^Kr + ^Be 
projectile fragmentation reaction [4^ . lYR of isobars are 
plotted in Fig. [21 which almost overlap except the lYR of 
/ =5. AB are calculated using the experimental binding 
energy of the nucleus in Ref. [i^l . The lYRs are fitted 
individually using Eq. ^ according to each /, and all 
the data are also fitted as a whole (shown as the solid 
line) . The fitted values of T are plotted in the inserted 
figure. Relatively similar T are obtained from / > 5 lYR. 
The line in the inserted figure represents T fitted from 
all the lYR data. 

To see T extracted from lYRs more systematically, 
lYRs in reactions of the lA GeV I24,i36xe (H, 790A 
MeV i29xe [IIL lA GeV "2,1243^^ [4^^ ^^^^ j4q^ 

"'^Ca and ^''Ni ^M,, are investigated. The results are plot- 
ted in Fig. [31 All the lYRs of isobars with different / are 
fitted as a whole and T obtained are given in each panel. 
lYRs in these reactions can be well fitted using Eq. ([5]). 
T and Afi obtained are plotted in Fig. [H The line in Fig. 
[HJa) is the average value of T in these reactions, which is 



28 • Relatively low T, which ranges from IMeV to 

29 3MeV, are found in these reactions. 

30 • The neutron richness of projectile affects T. T 

31 of a neutron-rich reaction system is higher than 

32 that of the neutron-proton symmetric reaction sys- 

33 terns when similar measurements are made, i.e., 
r(i3*^Xe) > T(i24xe), and T{^^^Sn) > T{^^^Sn). 

35 This is a similar phenomenon as the isospin depen- 

36 dence of the fragment yields measured in reactions 

37 of neutron-proton symmetric and neutron-rich re- 

38 actions [49j . The isotopic temperature (THeLi) was 

39 also found to increase when the projectile becomes 

40 more neutron rich using the isospin-dcpendent 

41 quantum molecular dynamics model [soj . 

42 • The mass of target affects T very slightly. T ob- 

43 tained from the ^®Ga -I- ^Be/^®^Ta reaction are very 

44 similar. A similar observation in T is made in the 

45 ^^Ni + ^Be/i^^Ta reaction. 

46 • The incident energy of the reaction, which ranges 

47 from 64A MeV to lA GeV, does not influence T 

48 very much. This occurs, as shown later using a 

49 modified statistical abrasion-ablation model anal- 

50 ysis, because T is essentially governed by the se- 

51 quential decay process. 
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FIG. 4. (Color online) (a) The fitted temperature (T), and 
(b) Afi from lYRs in reactions analyzed in Fig. [S] of 
tlie *®Kr and ^^^Xe reactions overlaps. The solid line in (a) 
represent the average value of T, and the line in (b) is the 
result of the linear fit between Afi and N/Z. 



side of Eq. ^ is rewritten as i?(A/i) = \nR{I + 2,I,A)- 
AS - Afi/T. In Fig. El i?(A/i) is plotted for typical re- 
actions as a function of AB for the isobar combinations 
of / = (7, 5) and / = (9, 7). T and A/i are taken from 
Fig. m For the ^''Ni reaction, R{Afi) values are slightly 
larger than those predicted by the average T. The mass 
range of the isotopes with / = (7, 5) and (9, 5) for ^"^Ni 
reaction is A =25 to 63. It was pointed out earlier that 
the analysis using A < 50 isobars has large uncertainty 
due to large Ato when T > 3MeV. Using the calcula- 
tion of Ato shown in Fig. [1] the experimental AB{Q) of 
isobars are temperature-corrected for T =-2 and 3MeV, 
which are labeled as AB{T)* . After the temperature 
correction, the correlation between the isobar combina- 
tions of / = (7, 5) and / = (9, 7) and the temperature- 
corrected AB{T)* are plotted in Fig. [S] It is easy to see 
that for the '^^Ca. reaction, T from AB{2)* and AB(3)* 
decrease to 1.56±0.12MeV and 1.27±0.10MeV, respec- 
tively. For other reactions, T from the temperature- 
corrected AB{T)* are only slightly modified. The cor- 
relation between the / = (7, 5) and I = (9, 7) lYRs and 
AB(3)* for the '^"Ni reaction overlap with AB/T. It can 
be concluded that T extracted from lYR of small mass 
is greatly affected by AB{T)* . 



T determined from lYRs of heavy fragments is lower 
than that from the Albergo isotopic temperatures of light 
fragments [TuHe) and slope temperatures (Tsiope) from 
the energy spectrum [l^ [^j , but is close to the isotopic 
temperatures involving heavier isotopes {TbcHc, Tubc, 
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BeLi: 



ThiLi |25[ and Tec ISll)- In a sequential decay 
process, light particles, such as p,n, and a, are emitted. 
When the light particle emission ceases, the fragment can 
still emit 7 rays and it further cools off. The extracted 
low temperatures associated with heavy fragments for the 
reactions in the wide incident energy range, indicate the 
dominance of the secondary decay effect on the temper- 
ature. One can also see the systematic trend that the 
temperatures from the lighter system tend to be higher 
than those of the heavier systems. Even though the MFM 
method provides a crude way to evaluate the tempera- 
tures, they are still useful in elucidating the entire process 
of the fragment production. 

, and A^ are associated with the properties of the 
emitting source but not the fragments themselves [l|, . 
Relatively large A^ were observed in this analysis for 
the reactions. A linear correlation is found between the 
correlation of A^u and N/Z of the projectile, which is 
shown as the solid line in Fig. HJ^b). The linear fitting 
result between Afi and N/Z reads y — (15.21 ± 5.82)2: + 
(-24.16 ±7.97). 

Rewriting Eq. ([S]), we get the following form. 



lni?(/ + 2, /, A)-AS- Afi/T = AB/T 



(6) 



with the left hand side involving lYR and A/i, which re- 1 
late to each reaction; the right hand side associates with 1 
AB/T of the isobars. For simplification, the left hand 1 
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FIG. 5. (Color online) The correlation between the i?(A/i) 
and the temperature-corrected AB{T)* of the I = (7,5) and 
I = (9, 7) lYRs. The solid symbols represent the results for 
the experimental AB{e); the crossed and open symbols rep- 
resent the results for the temperature-corrected AB{2)* and 
A_B(3)*, respectively. The lines are the fitting results using 
Eq. ([S]). The inserted figure shows the fitted values of T and 
A/i from the original lYR (ZnR(I-h2,I,A) - AS*). 



In Refs. [30, ImI, llS l52|, different temperatures 
(T =2.2, 6.0 and 9.5MeV) are used. Comparing to 
other works, a relatively low T =2.2MeV was used to 
estimate the binding energies of very neutron-rich cop- 
per isotopes [lOl in the 64A MeV *^Kr -|- ^Be reaction 
l43|. The staggering in the isotopic yield is minimized 
by introducing the approximation of back-shifted Fermi 



gas level density, or a parameter e. In Fig. [2l T and 31 
obtained from all the lYRs is T = 1.97 ± 0.10MeV32 
and —3.16 ± 0.26MeV, respectively. The dashed line 33 
represents T = 2.2MeV. In Ref. [s^, the equation 34 
used [Eq. (1)] to calculate the yield of a fragment 35 
(N, Z) is very similar to what we used in this article, 36 
Y{N, Z) = CA3/2exp[(Af/i„ + Z^J.p - F)/T], where F is 37 
the free energy, A/x = — 2.5MeV (/x„ = — 9.5MeV and 38 
^ip = -T.OMeV), T = 3/2 is different from the MFM of 39 
T = —3.6 ~ —2.2 (minus sign) in different / values [4l|. 40 
The very similar temperatures in Ref. [30| and this work 41 
indicate that temperature extracted from isobaric yield 42 
ratios is reasonable. 43 
Finally, we investigate the temperature of pre- fragment 44 
and final fragment in a modified statistical abrasion- 45 
ablation (SAA) model [H,!!!]- The SAA model can well 46 
reproduce the yield of fragments [i^ and was used to 47 
study the isospin phenomena in HICs [ssl - lssj . In brief, 48 
in the SAA model, the pre-fragment is calculated after 49 
the numbers of abraded protons and neutrons are known, 50 
which are determined by the nuclear-density distribution 51 
in the overlapping zone of projectile and target, and the 52 
nucleon-nucleon reaction cross section. Mean excitation 53 
energy of 13.3Ayl McV is assigned in the initial pre- 54 
fragments when AA numbers of protons and neutrons 55 
are removed from the projectile in the ablation-abrasion 56 
process. After the abrasion, the excited initial projectile 57 
nucleus undergoes the deexcitation process and forms the 58 
final fragment. The model description can be found in 59 
Rcfs. 49, 53, Sl. 
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FIG. 6. (Color online) (a) lYRs for the pre-fragments (crossed 
squares) and final fragments (solid squares) in the 140^4 MeV 
®*Ni -I- ^Be reaction of the SAA result, and those for the 
measured fragments (open circles), (b) T (circles) and Afi " 
(triangles) determined from lYRs of the prefragments (half- ™ 
full symbols), final fragments (full symbols) and the measured 
fragments (crossed symbols) in the 140A MeV """Ni -I- ^Be. (c) 8o 
and (d) are the same as that of (a) and (b), respectively, but 8i 
for the lA GeV "^Xe + 2°*Pb reaction. 82 

83 

The UOA MeV ^^^Ni 4- ^Be and lA GeV ^^sxe + ^ospb 84 
reactions are calculated. The lYRs of pre-fragments and 85 



final fragments in the Ni and -'^^^Xe reactions arc plot- 
ted in Fig. [nija) and (c), respectively. The lYRs of fi- 
nal fragments mostly overlap with those of the measured 
ones, while the lYRs of the pre-fragments have large dif- 
ference from those of final fragments and the measured 
ones. T and Afi determined from the lYRs with differ- 
ent / of pre-fragments and final fragments in the ^^Ni 
and ^^^Xe reactions are plotted in Figs. IHJb) andlBJd), 
respectively. T from lYRs of the final fragments are very 
similar to those from the experimental fragments, while 
T from lYRs of the pre-fragments are rather higher than 
those from the experimental fragments. It is shown that 
the drastic modification of T and chemical potential be- 
tween the pre-fragments and the final fragments results in 
the similar values in T and chemical potential extracted 
by the lYRs from the cold fragments in different reac- 
tions studied here. In other words, they are essentially 
governed by the sequential decay process. Therefore the 
yield of the cold fragments can be obtained by a sim- 
ple scaling of Eq. ([S]) for a variety of reaction systems 
in the wide incident energy range studied in this work. 
The extracted T plotted in Fig. Hl^a) show low T values of 
l~3McV. However, within that range they tend to corre- 
late with the projectile masses, that is, T decreases as the 
projectile mass increases. This may reflect the difference 
in T of the pre-fragments before the secondary decays. 
In the SAA analysis, the excitation energy of the pre- 
fragments are given by E* = 13.3AA MeV, where AA is 
the number of nuclcons removed from the projectile by 
the ablation-abrasion process. If the pre-fragment is in 
a thermal equilibrium, T will be given by T = ^JW^Ja 
and a = A/k, k is the inverse level density parameter. 
This leads to T = ^13.3/cAA/^. AA values are simi- 
lar for the reactions studied here. Therefore T becomes 
higher for the lighter pre-fragments which are produced 
more from the lighter projectile. However one should 
note that the extracted temperatures by the lYR from 
the cold fragments are significantly modified from that of 
the pre-fragments by the sequential decay process. The 
correlation in Fig. |4|a) merely reflects that the difference 
of the pre-fragments still sustains in some extent through 
the sequential decay process. 



IV. SUMMARY 

In summary, the temperature of a fragment after se- 
quential decay is studied using the isobaric yield ratio 
method in the framework of a modified Fisher model. 
The difference between the binding energy of the / > 5 
isobars at zero T [AB{Q)] is found to be valid for sub- 
stituting the value of AB{T) at low T. Relatively low 
T which range from 1 to 3MeV are obtained in different 
reactions. It is shown that T depends on the neutron- 
richness of the projectile. The mass of the target used 
affects T only slightly. The incident energy is found to 
affect T very little. A^ is found to depend linearly on the 
N /Z of projectile, i.e., larger A^ is found in reactions in- 
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86 duccd by more neutron-rich projectile. Due to the mass 22 

87 of the isobars of A < 50, an attempt was made to use 23 

1 the temperature-corrected AB{T)* . It is found that T of 24 

2 the ^^Ca reaction, in which isobars of A < 50 are dom- 25 

3 inant, is largely modified by the temperature-corrected 26 

4 AB{T)*, while T of other reactions, which involve iso- 27 

5 bars of larger masses, only arc slightly affected by the 28 

6 AB{T)*. 

7 The SAA model analysis for the UOA MeV ^''Ni-h^Be 

8 and lA GeV ^^^Xe+'^'^^Ph reactions revealed that the 

9 secondary decay process significantly modifies T and Ajj, ^' 

10 of the pre-fragments and governs those obtained from 

11 the cold fragments. This leads to similar lYR distribu- 30 

12 tions and temperature of cold fragment in reactions of 31 

13 different incident energies and different masses of projec- 32 

14 tiles. The SAA results also suggest that T from lYRs 33 

15 indeed reflects the actual physical temperature for the 34 

16 pre-fragments and the final fragments, although the lat- 35 

17 ter should be viewed as a sequenced temperature of the 36 

18 secondary decay process in conjunction with that of the 37 

19 primary process. Since the MFM method should be ap- 38 

20 plied to the initial hot nuclear matter in principle, the 39 

21 drawback of the application of this method to the fi- 40 



nal fragments is not directly to probe the characteris- 
tic nature of the initial stage, but mainly to probe that 
of the secondary decay processes. However, we believe 
that to elucidate the effect of the secondary process on 
the fragment yield is crucial to studying the nature of 
the primary emitting source, because all experimentally 
observed fragments have to go through this process. 
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